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ABSTRACT 

Several authors have reported that the dynamical masses of massive compact galaxies (M* > 
10 11 Mg, r e ~ 1 kpc), computed as Md yn = 5.0 o 2 r e /G, are lower than their stellar masses 
M*. In a previous study from our group, the discrepancy is interpreted as a breakdown of 
the assumption of homology that underlie the Md yn determinations. Here, we present new 
spectroscopy of six redshift z ~ 1.0 massive compact ellipticals from the Extended Groth 
Strip, obtained with the 10.4 m Gran Telescopio Canarias. We obtain velocity dispersions 
in the range 161-340 km 1 . As found by previous studies of massive compact galaxies, 
our velocity dispersions are lower than the virial expectation, and all of our galaxies show 
Md yn < M* (assuming a Salpeter initial mass function). Adding data from the literature, 
we build a sample covering a range of stellar masses and compactness in a narrow redshift 
range z « 1 . 0 . This allows us to exclude systematic effects on the data and evolutionary 
effects on the galaxy population, which could have affected previous studies. We confirm that 
mass discrepancy scales with galaxy compactness. We use the stellar mass plane (M*, a e , 
r e ) populated by our sample to constrain a generic evolution mechanism. We find that the 
simulations of the growth of massive ellipticals due to mergers agree with our constraints and 
discard the assumption of homology. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: evolution - galaxies: funda¬ 
mental parameters - galaxies: high-redshift - galaxies: kinematics and dynamics - galaxies: 
structure. 


1 INTRODUCTION 


The vast majority of early-type galaxies (ETGs) can safely be as¬ 
sumed to be in a steady stationary state. For these systems, the virial 
theorem holds, hence total potential and kinetic energies are bal¬ 
anced, such that 


GM ( v 2 ) t 
— ~ k ‘ — 


in 


where M is the dynamical mass, G is the gravitational constant, 
(r) T is the mean of the stellar distances to the centre of the system 
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at time r, (v 2 ) T is twice the kinetic energy of the system per unit 
mass at time r and liE ~ 2 is the virialization constant. 

Under the condition of structural and kinematical homology, 
we have (r) T oc r e and ( v 2 ) T oc o 2 \ we then expect that ETGs 
verify relationships such as M oc cr 2 r e and, wherever stellar mass 
dominates over dark matter, A/* oc a 2 r e ; here M* is the stel¬ 
lar mass (computed from population synthesis modelling), r e is 
the effective (projected half-light) radius and a e is the luminosity- 
weighted second moment of the line-of-sight velocity distribution 
(LOS VD; hereafter, velocity disp ersion). 

ICappellari et al.l (1 200(1. 1 20131) used Jeans and Schwarzschild 
modelling to provide a calibration of the dynamical mass for mas¬ 
sive elliptical and lenticular galaxies in the nearby Universe, such 
that 


© 2015 RAS 








Constraints on the evolution of massive galaxies 705 


Mdyn = 2 X M 1/2 = 5.0 —(2) 

Here Mdyn is defined as twice M \/ 2 , the mass enclosed in the iso¬ 
surface where half the light is emitted. Therefore, the ratio Md yn 
to M* provides a rigorous comparison of dynamical and stellar 
masses within the effective radius (because Md yn and M* are, re¬ 
spectively, twice the dynamical and stellar masses within the iso¬ 
surface where half the light is emitted). 

Dynamical masses of high-redshift galaxies, using equa¬ 
tion 0- have ofte n led to th e unphy sical resul t Md vn < M* 
(e.g. Stockton^hih&Marsgnl 201(J; Martine^MansojdaJ 20jl|; 


BarroetalJ 120141 : iHsu. Stockton & Shihl 120141 : IStockton et al.1 
20l4) . Significantly, the mass discrepancy is also found in nearby 


massive compact ellipticals jFerre-Mateu et al 

_ Addressing _ the _ mass dis crepanc y problem, 

IPeralta de Arriba et al.l d2014l , hereinafter |PdA14l) demonstrated 
that mass discrepancy does not depend on redshift, and instead 
scales with galaxy compactness, defined as the r e offset from the 
_stellar mass-size distribution of ETGs in the nearby Universe, 
derived an empirical scaling of Md yn with <r e and r e such 


that 


,, „ 3.6 0.35 

TF^dyn OC O e T e , 


(3) 


which implies a breakdown of ho molo gy. 

The galaxies analysed bv lPdAlj cover a wide redshift range 
(0 < z < 2.5) and come from diverse, heterogeneous sources. 
In this work, we analyse the Md yn -M* relationship over a narrow 
redshift range around 3 ~ 1.0. Hence, we exclude systematic ef¬ 
fects on the data and evolutionary effects on the galaxy population, 
which could plausibly have affected previous work. 

The core data set are 10 velocity dispersion measurements of 
9 extremely compact massive ellipticals at z ~ 1, taken with the 
10.4 m Gran Telescopio Canarias (GTC), 6 of which are new mea¬ 
surements that we describe in detail. We combine these data with 
velocity dispersions from the literature and draw a sample that cov¬ 
ers stellar masses 10 11 M 0 < M* < 10 12 M 0 and effective radii 
ranging from ‘normal’ to the most compact massive galaxies. This 
sample is ideal to show whether, at a narrow, cosmologically dis¬ 
tant epoch, the deviation between dynamical and stellar mass fol¬ 
lows the scaling relationship proposed bv lPdA 1 4j . We use a large 
z ~ 0 sample of massive galaxies to further test the predictions by 
comparing two snapshots in cosmic time. 

The paper is organized as follows. In Section [2] we describe 
the samples used in the paper: massive compact galaxies at z ~ 1 
(Section ITU the additional data at 2 ~ 1 (Section [2.21 and the 
nearby reference at z ~ 0 (Section 12.31 . In Section 12.1.11 we 
also report the recent M* values for our new six massive compact 
galaxies, extracted from the Rainbow Cosmological Sur veys data 
base (see lPerez-Gonzalez et al]l2008l : iBarro et al.ll2oTT3 lbh : while 
in Section ed.2l we describe the aperture correction for the velocity 
dispersions. In Section[3] we explain the spectroscopic observations 
(Section [3.1b . data reduction (Section [XT} and velocity dispersion 
measurements (Section lXX ). Accura te sky subtraction, done by im¬ 
plementing the method of iKelsonl j2003t ) within the REDUCEME 
package (Cardiel 1999), is demonstrated in Appendix 0 while 
the robustness of our velocity dispersion measurements is shown 
in Appendix [B] In Appendix [C] the reduced spectra and the fits 
are plotted. Appendix [D] gathers the structural parameters for the 


z ~ 1 sample. We compare dynamical and stellar masses in Sec¬ 
tion [4] confirming the dependence on compactness. In Section [5] 
we study the relationship between M*, r e and o e (the stellar mass 
plane), and use this information to constrain a generic evolution 
mechanism in Section [XT] In Section [5~2l we compare our results 
with numerical simulations of dry mergers. The discussion is pre¬ 
sented in Section [6] while the conclusions are given in Section [7] 
We adopt the concordance A cold dark matter cosmology (fl m =0.3, 
Ha=0.7, Ho = 70 km s -1 Mpc -1 ); at a = 1, 1 arcsec corre¬ 
sponds to 8.01 kpc. The stellar masses assume a Salpeter initial 
mass fu nction (IMF); where neces sary, we have used the relation¬ 
ships in lLonghetti & Saraccol ^20091) to convert stellar masses from 
the literature to this IMF. 


2 SAMPLES 

The core sample for this study is a set of nine massive com¬ 
pact galaxies at 3 ~ 1 for which our group obtained spec¬ 
troscopy on two runs with the Optical System for Imaging and low- 
Intermediate-Resolution Integrated Spectroscop}0 on the 10.4 m 
G TC. Data for four galax i es obt ained in the first run are described 
in iMartinez-Manso et al.l i201 ll . hereinafter MM11). The present 
paper fully describes the data from the second run, including sam¬ 
ple selection and characterization (Section [2.1b . spectroscopy and 
velocity dispersion determination (Section^- 

We expand this massive compact galaxy sample with data 
from the literature to create a z ~ 1 sample with a broad cover¬ 
age of the M*-r e plane (Section l2.2b . 

We also devise a reference z ~ 0 sample (Section [23l >. Galax¬ 
ies from both samples are similarly massive (M* > 10 11 M 0 ) 
and have spheroid-like surface brightness profiles (Sersic index 
n > 2.5). 


2.1 Massive compact galaxies at z ~ 1 

Target selection for both runs was similar. We selected pairs of 
galaxies from the catalogue of spheroid-like (n > 2.5), massive 
(M* > 10 11 M 0 ) galaxies in the Extended Groth Strip (EGS; 
iDavis et ai1l2007l) . built bv iTruiillo et alj d2007l . hereinafter |T07|) . 
This catalogue covers a redshift range from 0.2 to 2.0 and is com¬ 
plete in stellar mass down to 10 11 M 0 (the deta ils o f its stellar mass 
estim ations are presented in lBundv et alj|2006l and lConselice et aD 
120071) . Galaxies were selected in pairs given that we would use the 
long-slit mode of OSIRIS, and would follow the strategy of placing 
two targets in each slit. 

The selection criteria were: 

(i) To have a Sersic index n > 4.0. 

(ii) To be closer than 4 arcmin on the sky, to allow simultaneous 
observation with the long slit of OSIRIS. 

(iii) To have a spectroscopic redshift z ~ 1.0. 

(iv) To be near the compact end of the size distribution at each 
stellar mass. 

(v) To have apparent magnitude Jab < 22.2. 

Out of the sample, four galaxies were observed in the first run 
and they are described in MMllj. In the second run, we observed 


1 http://guaix.fis.ucm.es/~ncl/reduceme 


2 http://www.gtc.iac.es/instruments/osiris 


© 2015 RAS, MNRAS 453. 1704^7201 


































































706 L. Peralta de Arriba et al. 



07.60 s 07.40 s 14 a 19"’07.20 s 38.40 s 38.20 s 14 ,1 20"'38.00 s 06.60 s 06.40 s 14 , ’20"'06.20 s 


RA (J2000) RA (J2000) RA (J2000) 

Figure 1. Images of the six massive compact galaxies at z ~ 1, from the Hubble Legacy Archive. These images were taken with ACS through the F814W 
filter (/ band). They are combined images and were smoothed using a Gaussian filter with a standard deviation of 1 pixel. Each panel shows a galaxy at its 
centre and covers 5 arcsec on a side. Galaxy identifications are shown in the top left comer of each panel. A scale bar of 5 kpc is shown in each panel. Dotted 
lines represent the slits used in the spectroscopic observations. 


Table 1. Observed properties of our sa mple of six massive compact galax ies at redshift ~1. (1) Galaxy identifications from the DEEP- 
2 Galaxy Redshift Survey jDavis et ~al]|2003 1120071; I Newman et alj|2013l) . (2) Right ascensions. (3) Declinations. (4) Apparent I -band 
magnitudes in the AB system from lTruiillo et a lJ 1 20071 here inafte r I T07h . (5) Apparent A' s -band magnitudes in the AB system from47)7. 
(6) Effective (half-light) radii along the semimajor axis from lTO’t (7) Ellipticides from lTO't (8) Spectroscopic redshifts from the DEEP- 
2 Galaxy Redshift Survey (their absolute errors are <10 5 ). (9) Slit identifications of the spectroscopic observations. (10) Exposure 
times on target of the observations. (11) Rest-frame signal-to-noise ratios per angstrom in the spectra in the region of Ca II H and K (i.e. 
3900^1-000 A). 


ID 

RA 

Dec. 

I 

As 

fle 

e 

z 

Slit 

fexp 

(S/AOrest 


(J2000) 

(J2000) 

(AB mag) 

(AB mag) 

(arcsec) 




(h) 

(A- 1 ) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

13004276 

14 h 19 m 04?81 

+52°49'58764 

21.55 

19.21 

0.24 

0.42 

0.97533 

1 

6 

14 

12028173 

14 h 19 m 07?34 

+52°45'31 , /59 

22.06 

19.59 

0.20 

0.20 

0.95764 

1 

6 

13 

12024790 

14 h 18 m lT42 

+52°42'54789 

21.89 

20.00 

0.09 

0.43 

0.96572 

2 

4 

19 

12024422 

14 h 18 m 30!79 

+52 0 42'56'/65 

22.19 

19.82 

0.26 

0.44 

1.02874 

2 

4 

14 

13017337 

14 h 20 m 38?08 

+52°56'07'.'27 

21.20 

19.08 

0.31 

0.69 

0.97242 

3 

3 

15 

13018611 

14 h 20 m 06?37 

+52°55'33 , .'61 

22.15 

19.99 

0.17 

0.26 

1.07939 

3 

3 

12 


six galaxies. We chose to have one galaxy (ID 12024790) in com¬ 
mon with the sample of MM 11, to allow comparison of the velocity 
dispersions. 

In Fig. [I] we show Hubble Space Telescope !Advanced Cam¬ 
era for Surveys (ACS) images of the six galaxies of the second run. 
Their observed properties, as well as parameters for their spectro¬ 
scopic observations, are listed in Table[I] For the remaining of this 
paper, our sample will refer to the galaxies of the second run. 


2.1.1 Stellar masses 

For the galax ies of the second tun, we upgraded the stellar masses 
from the lT07l values to the data available in the Rainbow Cosmolog- 
ical Surveys data base (see Perez^Gonzalez et al.|2008j; Barro et alj 
1201 lal fbh. We chose the masses from Rainbow Cosmological Sur¬ 
veys data base because they have been computed usin g pho tom- 
etry with a broader spectral range (UV-FIR), while IT07I used 


only BRIJK. We used stellar mass es co mputed with spectro¬ 
scopic redshifts. These m asses assum e aSalgetM) 19551) IMF, mod¬ 
els from PEGASE 2.0 dFioc & Rocca-Volmerangel 1 9971) and a 
ICalzetti et alj d200 0f) extinction law. The upgraded values for these 
galaxies are on average 0.14 dex lower than the !t 07| values. This 
difference can be explained due to the usage of different models 
in their computation (see lMuzzin et alj|2009l) . Differences are less 
than lcr for each individual galaxy. 


These upgraded values for the galaxies of the second run have 
been included in Table [2] which summarizes derived properties of 
our sample of six compact galaxies. 

For the galaxies of the first run, IMMI ll report stellar masses 
assuming a Chabrier IMF; we have converted their values to a 
Salpeter IMF. After this conversion, the stellar mass determina¬ 
tions injMMll and Rainbow Cosmological Surveys data base are 
nearly but not completely identical: MM 11 assume a solar metal- 
licity, while the metallicity is not constrained in the stellar masses 
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Table 2. Derived properties of our sample of six massive compact galaxies at redshift ~1. (1) Galaxy iden¬ 
tifications from the DEEP-2 Galaxy Redshift Survey. (2) Sersic indices from ITO'A (3) Circ ular e ffective 
(half-light) radii com puted using the apparent effective radii along the semimajor axis (from lT07h . the el- 
lipticities (from lT07h and the redshifts (from the DEEP-2 Galaxy Redshift Survey); their uncertainties are 
10 per cent. (4) Stellar masses from |T07l ; their uncertainties are ^0.2 dex. (5) Stellar ma sses from Rain¬ 
bow Cosmological Surveys data base (see lPerez-Gonzalez et all2008tlBarro et al]|201 1 Jit) : in particular, 
the masses with the tag ZSPEC of this data base, which improve the fits using spectroscopic redshifts. (6) 
Statistical errors of the stellar masses from Rainbow Cosmological Surveys data base. (7) Mean velocity 
dispersions inside the slit aperture measured in 100 Monte Carlo realizations of each spectrum (errors 
are estimated with their standard deviations). (8) A per ture correction factor (computed following the pre¬ 
scriptions in the appendix B of Ivan de Sande et alJl2013h . (9) Velocity dispersions within the effective 
(half-light) radius. 


ID 

n 

r e 


M* 

AM* 

(7 

cr/(Te 

CTe 



(kpc) 

(10 11 M 0 ) 

(10 11 M 0 ) 

(dex) 

(km s -1 ) 


(km s -1 ) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

13004276 

5.19 

1.46 

3.5 

2.0 

0.09 

320 ± 15 

0.94 

340 ± 16 

12028173 

5.99 

1.42 

1.6 

1.2 

0.06 

214 ± 16 

0.94 

228 ± 17 

12024790 

4.47 

0.54 

1.2 

0.9 

0.03 

241 ± 14 

0.92 

261 ± 15 

12024422 

4.84 

1.57 

2.3 

1.2 

0.02 

225 ± 17 

0.94 

239 ± 18 

13017337 

4.33 

1.37 

3.6 

1.8 

0.5 

274 ± 23 

0.94 

291 ± 25 

13018611 

4.76 

1.19 

1.7 

3.0 

0.05 

151 ± 15 

0.94 

161 ± 16 



M* (10 u Me) 


z ~ 0 


SDSS NYU sample 


Z ~ 1 

0 This work 

■ Belli, Newman & Ellis (2014) 
Martinez-Manso et al. (2011) 
◄ Newman et al. (2010) 

van der Wei et al. (2008) 

& Blakeslee et al. (2006) 


Figure 2. The stellar mass-size distribution of the spheroid-like galaxies used in this work. The sources of the z ~ 1 sample a re given in the legend. The 
data point with a black edge represents the galaxy ID 12024790 also studied by I Martinez-Manso et all 1201 ll hereinafter lMMl ll) . The error bar cross on the 
bottom right corner represents the mean errors of our six compact galaxies at z ~ 1. The z ~ 0 sample is represented with a g rey-scale giving t he density 
of galaxies in the M*-r e plane. The dashed line shows the relationship that spheroid-like galaxies follow in the nearby Universe Ishen et al j2003l) . while its 
error bars show the dispersion of this relationship. 


from the Rainbow Cosmological Surveys data base. MM11| values 
are on average 0.11 dex lower than the values from Rainbow Cos¬ 
mological Surveys data base. Nevertheless, the differences fall well 
within the uncertainties of each determination. We will use the stel¬ 
lar mass values from MM11 for the galaxies of the first run in the 
remaining of the paper (except for the galaxy in common with the 
second run). 

In Fig. [2] we show the stellar mass-size distribution for the 
nine compact galaxies, together with that from the additional data 
described in Sections fc.Zl and ITTI Appendix[D]gathers the structural 


parameters for the z ~ 1 sample, including the stellar masses from 
the two runs of the core sample. 


2.1.2 Aperture correction of velocity dispersions 

In this paper, we will work with velocity dispersions integrated 
over an aperture of one effective radius. For the galaxies of the 
first run, MMlT| report velocity dispersions corrected to an aper¬ 
ture diameter of 1.19 kpc. We have taken their uncorrected values 
(around 6 per cent lower than the 1.19 kpc aperture-corrected val¬ 
ues) and performed a correction to one effective radius following 
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the prescriptions in the appendix B of Ivan de Sande et alj ( l2013h . 
These prescriptions take into account the effect of the point-spread 
function (PSF). To consider this fact is key due to the small ap¬ 
parent size of the galaxies at redshift 1. Also it is important to 
note that this correction can be used for apertures bigger than ef¬ 
fective radius; while this does not_occur for the aperture correc¬ 
tion eguations_pro£osed_by Ijorgensen. Franx & Kj aergaarjj[l99j) 


and 


Cappellari et alj d2006l) (see appendix B of lvan de Sande et al.l 


d2013l) for a detailed explanation). This correction increases around 
7 per cent their uncorrected values, i.e. the final values are similar 
to the 1.19 kpc aperture-corrected values. Appendix [D] gathers the 
structural parameters for the 2 ~ 1 sample, including the velocity 
dispersions from the first run of the core sample. 


2.2 Additional galaxies at z ~ 1 


We have obtained from the literature a sample of 18 massive ETGs 
at 2 ~ 1 (M* > 10 11 Mg, n > 2.5 and 0.9 < 2 < 1 . 1 ) with 
published effective radii, spectroscopic redshifts and velocity dis¬ 
persions. The sources of these data are: 


(i) iBelli. Newman & Ellisid2QI4l) : six galaxies. 

fii) iNewman et all j20 1 (t): one galaxy ( extracted from the com- 

pilatio nmade bvlv an de Sande et al. l20l 3jl _ _ 

(iii) Ivan der Wei et al] <1 20081) and iBlakeslee et al.l d2006h : 
11 galaxies (also extr acted from the compilation made by 
Ivan de Sande et akll2013l) . 


These publications offer the velocity dispersion values cor¬ 
rected to an aperture of one effec tive radius. They al s o foll ow the 
prescriptions in the appendix B of lvan de Sande et al.l <20131 ). 

The stellar mass-size distribution of these 2 ~ 1 data is shown 
in Fig. [2] All authors use our cosmology so no further corrections 
were needed. 


2.3 The nearby reference: massive galaxies from the SDSS 
NYU sample 

We have selected a sample of 53889 galaxies from the New York 
University Value-Added Galaxy Catalogue l lBlanton et al.ll2005h . 
which is based in the Sloan Digital Sky Survey (hereinafter SDSS; 
I York et aD EoOO). In what follows, we will refer this data set as 
SDSS NYU sample. The criteria for the selection of these galaxies 
were: (i) they have to be massive (10 11 Mg < M* < 10 12 Mg), 
(ii) they are spheroid like in > 2.5 in the r band), (iii) they 
have redshifts around the peak of the SDSS redshift distribution 
(0.05 < 2 < 0.11), (iv) they have reliable velocity dispersions 
(70 km s^ 1 < a < 420 km s -1 ) and (v) their physical sizes are 
reliable (0.3 kpc < r e < 30 k pc in the r ban d). 

Using equation 1 of ICappellari et all l l200dl . we have cor¬ 
rected the velocity dispersions from the fixed aperture (3 arcsec) 
of the SDSS fibres to the expected within one effective radius. 

The stellar mass-size distribution of these data is shown in 

Fig. 0 


3 DATA 

In this section, we explain the spectroscopic observations, data re¬ 
duction and velocity dispersion measurements for the second run 
of six massive compact galaxies described in Section fTTl 


3.1 Spectroscopic observations 

Long-slit spectra for the six galaxies were obtained using the 
OSIRIS instrument at the 10.4 m GTC telescope. Observations 
were carried out in queue-scheduled service mode in 10 nights be¬ 
tween 2012 April and August. We used the grating R2500I, which 
covers wavelengths from 7330 to 10000 A, corresponding rest- 
frame wavelengths at 2 = 1.0 are 3665 A < A re st < 5000 A. 
The slit width was 0.8 arcsec, which delivered a resolving power 
of R = 2050 at A = 8000 A as measured from the width of 
the arc lines, which translates into an instrumental resolution of 
<j mat ~ 64 km s -1 . We binned the CCD 2x2, yielding a sampling 
of 0.254 arcsec and 1.36 A in the spatial and spectral directions, re¬ 
spectively. Hence, we sampled the spectral resolution element with 
2.9 pixels. Observations were carried out in dark time with seeing 
<0.9 arcsec FWHM (full width at half-maximum). 

Observations were broken in observing blocks of 1 h, each 
containing two 30 min exposures dithered by 7 arcsec along the slit. 
Total on-source integration for the three slits was 16.5 h, but 3.5 h 
yielded low counts, probably due to high clouds, and were unused. 
Exposure times on target for slits 1, 2 and 3 were 6, 4 and 3 h, 
respectively. Spectrophotometric standards from the Isaac Newton 
Group of Telescopes librar}@ were observed each night. 


3.2 Data reduction 

We reduced our data following the standard steps (bias subtraction, 
flat-field correction, sky subtraction, cosmic ray removal, wave¬ 
length calibration, C-distortion correction, extinction correction, S- 
distortion correction, extraction and relative-flux calibration) using 
REDUCEME jCardieli 1999t> . a reduction package for long-slit spec¬ 
tra. Sky subtraction in spectra redwards of 7000 A is particularly 
difficult because the profiles of strong OH telluric lines, sampled 
by a few pixels only, get affected by spectral rectification yielding 
strong sky-subtraction residuals. We implemented in REDUCEME 
the algorithm proposed bv lKelsoiil ( 120030 : taking advantage of the 
gradual shift of the wavelength solution along the slit, we build a 
grossly oversampled sky spectrum, which we then rebin using ap¬ 
propriate shifts for each CCD row. In Appendix [A] we describe the 
software tools developed to use this method and show examples 
of its effectiveness. Source extraction was carried out via straight 
sum within an aperture of 9 pixels (2.286 arcsec). Signal-to-noise 
ratios measured in the 3900^1000 A rest-frame wavelength range 
from 12 to 19 per angstrom; they are listed in column 11 of TableQ] 
Relative-flux calibration and correction for telluric absorption was 
carried out by using the spectrophotometric standard stars. 

The reduced spectra of the six galaxies are plotted in Ap¬ 
pendix 0 together with the spectral fits described in Section 1331 
The spectra cover the rest-frame range 3800-4500 A. Call H 
and K, the G band and several Balmer lines feature prominently, 
with varying relative intensities, indicating a mixture of old and 
intermediate-age stars. In addition. Mg I 3829, 3838 A is visible in 
some of the spec tra. The type of spectral features is similar to that 
found by MM11. Given the blue cut-off of our spectra, the presence 
of the forbidden emission line [O II] 3727 A can only be checked 
on two of our galaxies; no firm evidence is found on any of the two 
galaxies (IDs 12024422 and 13018611). 


3 http://catserver.ing.iac.es/landscape/tn065-100/workflux.php 
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M* (10 11 Mq) M* (10 11 M 0 ) 


Figure 3. Stellar versus dynamical masses for the two redshift samples. Dashed line shows the identity relationship (AT* = A/ ( | yn >■ Dotted lines below/above 
the dashed line indicate how many times (in powers of 2) the stellar mass is greater/smaller than dynamical mass (i.e. M* = 2 Z A/,i yn where i is an integer). 
Left-hand panel shows the z ~ 0 sample. For this sample, the data have been binned depending on their dynamical and stellar masses. The size of each 
symbol in the left-hand panel scales with the number of galaxies in each bin. For clarity, bins with fewer than 10 galaxies have been omitted. In the right- 
hand panel, which contains galaxies with z ~ 1, symbol shapes as in Fig. [2] In both panels, the colour of each symbol represents the compactness indicator 
r e/ r Shen(JVf*)> d e fi ne( l as the r e offset from the stellar mass-size distribution of ETGs in the nearby Universe. The error bar cross on the top of the right-hand 
panel represents the mean errors of our six compact galaxies at z ~ 1. 


3.3 Velocity dispersion measurements 

We computed velocity dispersions using the penalized pixel fitting 
(pPXF) method of Icappellari & Emselleml (120041) . This method 
characterizes the LOSVD by decomposing a spectrum as a lin¬ 
ear superposition of templates which have been shifted consid- 
ering the Doppler effect. We followed the steps described in 
iFalcon-Barroso et all d201 ill , fitting the LOSVD with the first two 
Gauss-Flermite moments and including the computation of uncer¬ 
tainties using Monte Carlo realizations. 

Our spectral template library comprised simple stell ar pop¬ 
ulations (SSPs) from the PEGASE-HR library ( me Borgne et ali 
1 20041) . This library was selected because its high spectral resolu¬ 
tion (0.55 A FWHM) allows us to take advantage of the quality of 
our spectra (which have a resolution of ~2 A FWF1M at rest frame). 
We restricted the choice of templates to SSP ages younger than age 
of the Universe at 2 = 1. 

The use of SSPs was selected by the fact that SSPs have a 
more physically motivated combination of line strengths than com¬ 
binations of stellar templates obtained by minimizing \ 2 ■ This is 
important given that Call H and K can bias velocity dispersion 
determinations because their intrinsic width varies with stellar at¬ 
mospheric temperature dGebhardt et alj2Q03h . Because our spectra 
(Appendix [C]l show prominent Balmer lines in addition to Call FI 
and K, the use of SSPs prevents the fitting programme against giv¬ 
ing a solution that depends strongly on the relative weights of the 
stellar template mix. We show in Appendix [B] that the spectral fits 
are particularly sensitive to the choice of stellar templates. 

The final fits are shown in Appendix [C] Although the faintest 
targets have higher relative sky-subtraction residuals, the main 
spectral features can be detected by visual inspection, and properly 
identified by the fitting programme. 

The velocity dispersions are listed in Table [2] together with 
derived properties of our sample of six massive compact galax¬ 
ies. The robustness of these results is explained in Appendix m 


Table [2] also provides the velocity dispersions corrected to one 
r e apertures. This correction follows the prescriptions in ap¬ 
pendix B of Ivan de Sa nde et af] l|2013l), whic h ad apts the cor- 
rections propose d b yjjorgensen, Franx & Kiaergaardl (119951) and 
ICappellari et ali HoOfj) to higher apertures and takes into account 
the PSF effects. 


4 DYNAMICAL VERSUS STELLAR MASS: 

DISCREPANCY GROWS WITH COMPACTNESS 

In this section, we compare the dyn amical mass computed as 
ATdyn = Ka^re/G with K = 5.0 i lCaPDellari et alj|2006l) with 
the stellar masses derived using stellar population techniques. 

In Fig. [3] we plot dynamical mass versus the stellar mass for 
our two redshift samples. We extract two results from this figure. 
First, in both panels we find galaxies below the dashed line. This 
indicates that the unphysical result AT* > ATdyn occurs not only at 
2 ~ 1 but also at 2 ~ 0. 

Our second result relates to the dependence on galaxy com¬ 
pactness. The colour of each symbol in Fig.[3]has been used to in¬ 
dicate the compactness indicator defined as re/)"shen(AT*), where 
t*shen (AT*) is th e mass-si ze relationship for galaxies in the nearby 
Universe from lshen et al.l ( t2003h , written as 

r Shen(AT*) = r Shen ^ 1Q11 ^ j (4) 

where rg hen = 3.185 kpc and p = 0.56. T hese values come from 
adapting the relationship in Shen_etjd](|2003b to a Salpe ter IMF fol¬ 
lowing the equations from Longhetti & Saracccl l2009t) . The com¬ 
pactness indicator measures the r e offset from the stellar mass-size 
distribution of ETGs in the nearby Universe. 

In the left-hand panel of Fig. [3] which shows the 2 ~ 
0 galaxies from the SDSS NYU sample, we can see that red 
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colours (low r e /rshen(Af*)) populate the area below the dashed 
line. This indicates that mass discrepancy in the nearby Universe 
grows with compa ctness. Observational errors in r e (15 per cent; 
ICebrian & TruiillJ 120141) translate into a 15 per cent scatter in 
Mdyn at a constant M*, which is only one third of the observed 
scatter of Md yn at a constant M*. This corroborates that the trend 
of low Mdyn with compactness in Fig.[3]is real and not a result of 
error propagation. 

At z ~ 1 (right-hand panel of Fig©, again, the distribution of 
the symbol colours follows a pattern such that more compact galax¬ 
ies show stronger mass discrepancies. It is worth noting that in this 
redshift we find higher discrepant values due to the strong size evo¬ 
lution of this type of galaxies (e.g . Daddi et al. 2005; Trujillo et al. 


lution ot this type ot galaxies (e.g. Uaadi et at. 2UU3; iruullo et al. 

20061; Longhetti et al J2007I: Toft et alj2007 : T07 : Zirm et al J2007I : 


Cimatti et al]|2008l;lBuitrago et al.ll2008h. 

Hence, our two samples covering narrow redshift ranges, and, 
therefore, largely free from redshift-dependent syste matics in the 
mass determinations, confirm the result from |PdA14 . namely, that 
mass discrepancy grows with galaxy compactness. 

We emphasize that, had we chosen a different IMF such as 
Chabrier, we would have alleviated the problem of galaxies show¬ 
ing Mdyn < M*, but would not have solved it_The range of 
IMF variations propos ed by authors such as iTreu et alj hold ). 
Cappellari et al. (2012) or lFerreras et alJ d201 3h would reduce our 
stellar masses by at most a factor of 2. We would still have sev¬ 
eral galaxies at the low-redshift panel showing the mass anomaly. 
However, at high redshift, selecting other IMF will be insuffi¬ 
cient to solve the problem. Furthermore, the authors mentioned 
above associate bottom-heavy IMFs with massive (high-velocity- 
dispersion) galaxies, which supports the choice of a Salpeter IMF 
made in this paper. This ex pectation agrees with the results from 
iMartln-Navarro et al J d20 15h . which found a bottom-heavy IMF for 
massive ETGs at z ~ 1. 


4.1 Weak homology: non-homology from differences in light 
profiles 

Several authors have proposed that a deviation from a constant 
value for K could be originated by the different light profiles that 
galaxies show. Using one-component , sph erical, non-rotating and 
isotropic models, iBertin. Ciotti & Del Principel d2002l) derived an 
equation where K dep ends on Sersic index n (from 6.6 to 2.2 in the 
interval n t (2.5, 8)1. ICappellari et aD d2006l ) found also a similar 
expression assuming spherical isotropic models with Sersic profile 
(from 6.9 to 3.8 when n varies from 2.5 to 8). Nevertheless, these 
authors fou nd a quite c onstan t value for their observed galaxies 
(5.0±0.1). lfavlor et alJd2010l) . comparing also stellar masses from 
stellar population techniques with virial mass estima tors, reporte d 
that fo llowing the prescriptions of IBertin. Ciotti & Del Principg 
( 12OO2I) improved the correspondence between stellar and dynam¬ 
ical masses. 

Could this effect remove the trend between M*/Md yn and 
t'e/?’shen(M*) reported in Section [4]? Fig. [4] addresses this ques¬ 
tion showing the relation between and r e /r'shen(M*), 

being MJyn the dynamical mass computed with a K which fol¬ 
lows the dependence on n modelled by lCappellari et al1 ( i2006l ). i.e. 

K(n) = 8.87 - 0.831n + 0.0241n 2 . (5) 

Thanks to this figure, we can check that the trend between 
M*/MJ^® k and r e /rshen(M*) rem ains. The trend also remains 
if we use the expression proposed bv lBertin, Ciotti & Del Principel 

d2002l) . 



Figure 4. Correlation between and the compactness indica¬ 
tor r e /rghen(M*). being the dynamical mass computed with the 

dependence of K on Sersic index n given by equation 0 . Symbols for 
the 2 ~ 1 sample as in Fig. [2] The solid black line shows the data from 
the 2 ~ 0 sample after grouping them in their 1000-quantiles of their 
r e/rghen(M*) distribution; the width of this line shows the error of the 
mean in each 1000-quantile. The grey area shows the unphysical region 
where M+ > Mdyn■ The dash-dotted line corresponds to a stellar-to- 
dynamical mass ratio_of 0.7, derived from lensing of nearby massive el¬ 
lipticals bv lGavazzi et all 12007) . 


5 STELLAR MASS PLANE 


In this section, we study the relationship between M*, r e and cr e . 
We fit a plane of the form M* oc cr“rg, with a and b real numbers, 
to the z ~ 1 sample. The hypotheses of virial equilibrium and 
homology imply a plane of the form 


M* = 


where 


K 


2 

OVe 


K is a constant 


(6) 


(7) 


and fdyn = M*/M d yn may also be considered as constant, due 
to a random or weak dependence on M», r e o r a e . In the fol¬ 
lowing, we assume /d yn = 0.7 iGavazzi etah 2007 f This as¬ 
sumption agrees with the results reported by Tavlor et all d2010h . 
which found that /d yn depends very weakly on stellar mass at fixed 
Sersic index n and varies little for the range of n of our sample 
(~0.1 dex). As noted in Section |T] several works have found that 
equation © applies t o nearby, massive ellipticals and lenticulars 
(e.g. lCappellari et al j20061). Because equation © leads to unphys¬ 
ical dynamical masses, PdA14l proposed a modification of equa¬ 
tion ©. where K is allowed to vary with the compactness indicator 
fe/rshen(AT*), such that 


K = 


fShen (M* 


( 8 ) 


where k is a constant. It is worth noting that this parametriza- 
tion allows also the variation of /d yn with the compactness indi¬ 
cator fe/rshen(M*), but including it as part of the K variation. 
These authors choose this option because they argued that the de¬ 
viation from homology cannot be produced mainly by the variation 
of/ d * yn (cf.fig.9oflMA14l). 

Our six compact galaxies add essential data to test for vari¬ 
ations of K with compactness. In Fig. [5] we show the K value 
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r e /r e (M*) 


Z 1 sample of this work: 
a = -0.63 ±0.17, k = 6.6 ±0.7 
Peralta de Arriba et al. (2014): 
a = -0.81 ±0.05, k = 6.0±0.2 
Homology with the K value 
found by Cappellari et al. (2006): 
a = 0 with K = k = 5.0 


Figure 5. K value versus the compactness indicator r e /rc^ len (M,) (i.e. the distance to nearby stellar mass-size scaling relationship). Sym bols a s in Fig. [2] 
Red solid and green dotted lines show the prediction of virial theorem with homology (considering the value K = 5.0 found by lCappellari et al]|2006l) . and 
the prediction of lPdAl j . The blue dashed line shows the fit of t he z ~ 1 sample. Black arrows show expected evolution up to z = 0 of the six massive 
compact galaxies if they would follow the equations proposed by[Tapia et alJ 120id see Section [5©. The error bar cross on the bottom left comer represents 
the mean errors of our six compact galaxies at z ~ 1. 


versus the compactness indicator fe/rshen(M*). The red solid line 
represents the predictions of virial theorem with homology (equa - 
tion© using the value K = 5.0 found bylCa ppellari et ah ibood). 
while we used the green dotted line to plot the predictions of |PdA14l 
(equation©, which proposed a = —0.81±0.05 and k = 6.0±0.2. 
All of our mass_ive compact galaxies except one follow the pre¬ 
dictions ol lPdA 141 while one of them shows a higher value of K 
than both predictions. Its extreme K value is caused by a low a e 
(161 ± 16 km s -1 ); K drops slightly (K = 40.3) but remains well 
above the PdA14 relation, if we use the velocity dispersion from 


iFemandez Lorenzo et al J J20I iL a e = 196 ± 6 km s *). 

Fitting equation © to all of our objects at z ~ 1 , we obtain 


a = -0.63 ±0.17 


(9) 


and 


where Ml, rl and al denote the initial stage of an individual galaxy 
before the interaction of the mechanism, and the superindex ‘f’ is 
analogously used to indicate the final stage. We choose such scal¬ 
ing relations because of their simplicity and because we expect the 
changes of r e and <r e to be strongly affected by the changes in mass. 

The geometric interpretation of the above equations is that the 
mechanism driving size evolution moves galaxies in a well-defined 
direction in the three-dimensional space (log M*, log a e , log r e ). 

If we consider that this movement has to happen inside a plane 
of the type defined by equations © and ©, we obtain a relation¬ 
ship between the exponents: 


p + 2E - 1 

v-P 


(13) 


k = 6.6 ± 0.7, (10) 

where the errors have been computed using the bootstrap re¬ 
sampling technique. This is lcr compatible with the results from 
|PdAl4 but it is inconsistent with the prediction of equation © (a 
is 4cr from 0). This fit has been represented in Fig. [5] with a blue 
dashed line. Note that in Fig. ©he regression to 2 ~ 1 (equations[9] 
andl 1 0b matches the z ~ 0 sample despite the latter being excluded 
from the fit. 


5.1 Evolutionary constraints of stellar mass plane 


In this section, we study the constraints that the stellar mass plane 
implies for the mechanisms which drive size evolution in this type 
of galaxies. 

We focus on a generic mechanism in which the changes of r e 
and a e scale as power laws of the mass growth: 



(ID 


( 12 ) 


where 77 is the slope of the mass-size relation for early types in the 
nearby Universe (equation©. Under the virial theorem and homol¬ 
ogy (i.e. a = 0 ) a simpler relationship holds: 

£ = (14) 

Additionally, relaxing the hypothesis of homology to allow a de¬ 
pendence of /j on M* (i.e. oc M*), it leads to a shifted 
analogous relationship: 


In the following, when we consider this last case, we will assume 
the value S ~ —0.2 (which comes from parametrizing the variation 
of ff n reported bv lTavlor et ali J20ld) with M*). 

In Fig.© we show the restrictions on the p-E space that sev¬ 
eral stellar mass planes imply: (i) the plane predicted by virial theo¬ 
rem and homology (red solid line), (ii) the plane predicted by virial 
theorem and a relaxed homology where /±„ oc A /* -0 ' 2 (red dash- 
dotted line), (iii) the plane proposed bv |PdA14l (green dotted line) 
and (iv) the plane obtained fitting the 2 ~ 1 sample of this work. 


© 2015 RAS, MNRAS 453. I704H7201 












































712 L. Peralta de Arriba et al. 



P 


DRY-MERGING SIMULATIONS 


i-fn Tapia et al. (2015) 



Figure 6. Stellar mass plane constraints on exponents p and S for a generic mechanism of evolution of the type r\jr\ = (Ml/M+J, <j\ = 

(M^/M^) . Blue dashed and green dotted lines show the constraints from a exponents obtained from the stellar mass plane proposed in Section [7] of 
this work and in 11'dA I -1 . The shadow regions around these lines cover the ler errors of these exponents. The red solid line shows the constraints from virial 
theorem and homology prediction. The red dash-dotted line shows the constraints from virial theorem and a relaxed homology where stellar-to-dynamical 
mass ratio depends on stellar mass (/£ °c Black data cross is the result found in the simulations of lTapia et all 120151) . The results from the 

simulations of Hilzet al. I20i2t > are represented as indicated in the legend: 1:1/1:10 label corresponds to simulations of major/minor mergers, the label 'DM’ 
indicates the galaxies are embedded in matter haloes (or not if it is crossed out), the label 'AM' refers to an orbit with angular momentum (or without angular 
momentum if it is crossed out), and the presence of the suffix ‘-c’ means that the merging satellites were compact (i.e. satellites sizes were 50 per cent of the 
sizes assumed in the other simulations). 


5.2 Numerical simulations of dry mergers 

For massive compact galaxies, the cosmological evolution on the 
(AT*, cr e , r e ) space due to their merger history can be extracted from 
numerical merger simulations. ITapia et alj <201511 extracted merger 
histories from cosmological hydrodynamical simulations and re¬ 
produced them from a = 2.5 to 2 = Oat high resolution, using 
TV-body techniques. Initial primary galaxies had M* ~ 10 11 Mq 
and r e = 1 kpc, while the merging secondaries had radii following 
the observational mass-size relationship at the epoch of the merger. 
Their simulations show that typically, a massive compact galaxy 
grows by a factor ~2 in mass and a factor ~4 in r e from 2 = 2.5 
to 2 = 0. 

The changes of r e and a e with mass growth measured in 
the simulations, fitted with equations GB and 03. yield p = 
1.94 ± 0.09 and E = 0.06 ± 0.02. These values are plotted in 
Fig. [6] There is a close agreement of the_results of merger simula¬ 
tions with the stellar mass planes of lPdAljj and that found in this 
work (equations [9l and ITOt . In contrast, the results of the merger 
simulations fall far from the prediction of the virial theorem and 
homology (even when this last assumption is relaxed to allow vari¬ 
ations of /E as function of M*). 

iHilz et alJ i 20121) simulated 10 series of mergers with mass 
ratios of 1:1 and 1:10, for a range of galaxy models (with or with¬ 
out dark matter haloes; varying the internal densities of the sec¬ 
ondaries) and galaxy orbits (with or without orbital angular mo¬ 
mentum). We have modelled the effects of these simulations with 
equations 03 and 03- Thi s allo ws to include the results of 10 
simulations from lHilz et alj 13)121 ) in Fig. [6] This figure shows that 
the action of major or minor mergers gives movements in the space 
(M*, (T e , r e ) which disagree with the homologous virial predic¬ 
tions, although major mergers and one of the minor-merger cases 
(mass ratio 1:10 without dark matter and with compact satellites 
which orbits do not have angular momentum) could be fitted with 
a deviation from homology where f£ yn cx: AT*. We can check also 


that most of the simulations (8 out of 10) fall within the ler error of 
the stellar mass plane proposed in this work. 

The comparisons with /V-body mergers described here show 
that the evolution of r e and a e as AT* grows due to mergers moves 
massive galaxies out of the homology plane, along the mass plane 
proposed here. Nevertheless, a deviation from homology varying 
stellar-to-dynamical mass ratio with stellar mass can match the 
major-merger simulations and a particular case of minor mergers. 


6 DISCUSSION 

It has been amply recognized that velocity dispersions contain 
key information to constrain the evolutionary processes linking 
massive compact galaxies at high redshift with nearby massive 
ellipticals: a simple transformation between equilibrium states 
by injection of dynamical energy predicts velocity dispersions 
around 500 km s _1 for galaxy masses ~ 10 11 Mq, while galaxy 
growth via mergers p redicts velocity dispersions to hardly evolve 
I lHopkins et alJEoioi ). Our measurements add to the growing body 
of data that shows velocity dispersions (<r e in the range 200- 
350 kms -1 ) for massive compact ellipticals. To our knowledge, 
there remains one single case of a massive passive galaxy with 
reported velocity dispersion above 500 kms -1 , namely, a 2 = 
2.186 gal axy wi th = 3 .2 x 1 0 11 Mq and r e = 0.78 kpc, 
for which Ivan Dokkum. Kriek & Franxl J2009I) reported a value of 
5102gf kms -1 . 

The comparison of our velocity dispersion measurements 
with those from other authors (see Appendix 0 indicates that 
uncertain ties in jvelocity dispersions could be underestimated. 
Ide Bruvne et al .1 120031 ) advised that the usage of Monte Carlo sim¬ 
ulations with white noise for estimating the errors can lead to values 
which are almost a factor of 2 lower than realistic errors. Although 
we have tried to improve our errors rejecting the strong sky resid- 
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uals in the fits and in the amplitude of the white noise for Monte 
Carlo simulations, the dependence between pixels introduced in the 
reduction steps should be the dominant factor for the underestima¬ 
tion of the errors. 

The explanation of the tilt of Fundamental Plane is a cl as- 
sic debate of extragalactic astronomy (e.g. lTruiillo. Burkert & Bell 
l200i: R enzinill2006l) . Two options (or a combination of both) are 
possible: stellar population effects, or deviation from homology 
(where we include also variations of dark matter fractions). The 
stellar mass plane used in this work has the advantage of replacing 
the luminosity by the stellar mass (from stellar population tech¬ 
niques). This allows to address the problem in a space where the 
three variables are dynamically connected. Therefore, the existence 
of a tilt in this space implies the violation of homology. 

The comparison with numerical simulations of dry mergers 
has worried us about some interpretation issues. For example, in 
the simulations from Hilz et al] d2012l) , although the break of ho¬ 
mology happens in our observables (M*, cr e , r e ), it does not occur 
in a ‘theoretical’ space where the variables are total bound mass, 
mean square speed and gravitational radius. This can be seen as an 
obviousness, because it only means that the equilibrium of ener¬ 
gies predicted by virial theorem is satisfied. Flowever, this advises 
that the assumption of homology is not possible for translating dry- 
merger simulation results from theoretical variables to observables 
(or even between different observables). 

We have also compared our results with numeri¬ 
cal _simulations_ of puffing-up model. The simulations of 
iRagone-Figueroa & Granatol OOl li) predict that the three- 
dimensional stellar velocity dispersion and the half-mass radius 
are related by 0-30 oc r^ 2 , with j 3 varying from —0.49 to —0.83 
depending on the mass loss from gas ejection and the ejection time 
chosen in the simulation. As fi is not always —1/2, this indicates 
that homology can also be violated in puffing-up simulations. 
Our stellar mass plane predicts a law a e oc r ~ 019±0 09 a t fixed 
stellar mass, while the plane from |PdA14l gives <r e oc r~° 10±0 03 . 
Assuming an ‘ad hoc homology’, i.e. r e oc rq/2 and cr e oc 03D, we 
could conclude that our stellar mass plane rejects a big contribution 
of puffing-up to the size evolution of ETGs. Nevertheless, when 
the break of homology is detected, this type of assumptions can 
lead to wrong results. Therefore, a puffing-up simulation which 
results where given in terms of observables would be useful to 
reject definitely the puffing-up mechanism. 

Fig. [6] shows the importance of constraining observationally 
the stellar mass plane. The predictions from minor-merger evolu¬ 
tion (Hilz et al.! 20121) fall well within the stellar mass plane pre¬ 
sented in this paper, while they deviate from the mass plane derived 
in lPdA14l While the set of merger models with which we compare 
are quite restricted, Fig. [6] shows that this type of analysis has the 
potential of allowing to distinguish between major and minor merg¬ 
ers as the dominant channel for the mass and size growth histories 
of today’s giant galaxies. While minor mergers havejtften been 
claimed for the growth of massive comp act galaxies (Honki ns~et al] 
120091 : 1 Hilz et alJEoIl iTania et alj|200) . a growth based on 1:1 
to 1:5 mergers is in fact suggested by observations of the satel- 
lite distribution around massive galaxies at different re dshifts (e.g. 

iFerreras et al.ll2014l : iRuiz. Trujillo & Marmol-Oueraltdl20f4h . 


7 CONCLUSIONS 

The main conclusions of this paper are: 

(i) Our measurements of velocity dispersions for six massive 


compact galaxies do not reveal the high values expected from a 

law M* OC CTefe. 


(ii) The discrepancy between stellar and dynamical masses com- 
putedjisMdyn = Ka^re/G with K = 5.0 follows the predictions 
of |PdAl4 i.e. it scales with galaxy compactness. This implies a 
breakdown of homology. 

(iii) For our sample of massive galaxies at 2 ~ 1 which in¬ 
cludes our six compact galaxies, we find a relationship between 
Mi,, errand r e compatible with the alternative scaling law proposed 
bv lPdAlj . This result is compatible with a snapshot at 2 ~ 0 with 
galaxies from SDSS. 

(iv) The relationship between M*, a e and r e is compatible with 
numerical studies of the growth of massive ellipticals due to a mix¬ 
ture of minor and major mergers that have a cosm ological frame¬ 
work like the one realized byiTaniaetaOj 2(M5j). The numerical 
simulations of dry mergers from Hilz et a l. d2012h also predict a 
break of homology in the same direction that our stellar mass plane 
constrains. 
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APPENDIX A: SKY SUBTRACTION IN THE DATA 
REDUCTION 

The most critical step of our data reduction was the sky subtraction. 
The reason of this trouble is based on the presence of many narrow 
lines of sky emission on the near-infrared. In order to try to get the 
best possible results on this step, we developed the new software 
tool PYKELSAME, which implements the sky-subtraction method 
proposed bv lKelsonl J2003h in the environment of REDUCEME. 

In addition, in the development of PYKELSAME we found the 
necessity of coding other two auxiliary tools: REDUClO and PYS- 
DISTOR. REDUClO is a simple module which allows to read and 
write files with the format of REDUCEME in a PYTHON environ¬ 
ment (REDUCEME has its own format because it was designed in 
the 90s and then the speed of reading and writing the data played 
an important role). The function of PYSDISTOR is to trace the spec¬ 
trum of faint targets (or any) in the CCD. The algorithm of this last 
tool is simple: to trace the spectrum adding several channels (i.e. 
pixels along the spectral direction) to make easier to detect a faint 
object. 

Once we developed these two auxil iary tool s, we implemented 
the sky-subtraction method proposed bv lKelsonl I I2003D . Fig. lAll il- 
lustrates the idea of this method: once we have characterized the C- 
distortion of the CCD (but not corrected it), we can use all the rows 
of the CCD to build a high-resolution spectrum of the sky. Fig. [ATI 
shows how by using every pixel of the image without rebinning 
(left-hand panel) we can build a high-resolution spectrum (black 
thin line on the right-hand panel) compared with the spectrum ex¬ 
tracted from a single row (red thick line on the right-hand panel). 
It is worth noting the presence a cosmic ray in the high-resolution 
spectrum of this figure: it is an example of an issue which has to be 
solved in the construction of the sky model. 

Fig. IA2I illustrates the process of the construction of the sky 
model. In top panel, we show the first step: the CCD region to be 
considered in the modelling. We checked that our results improve 
selecting a CCD region around the target but excluding the region 
affected by the target spectrum (we detected the target spectrum 
using PYSDISTOR); the dashed lines in the top panel show the ex¬ 
cluded region. The middle panel shows the high-resolution spec¬ 
trum built from the CCD region using the technique explained in 
the previous paragraph. It is worth noting the existence of cosmic 
rays in this spectrum. To reject these points, the sky model is con¬ 
structed as follows: we clean the high-resolution spectrum using 
a median filter (or any other percentile if it was desired), and we 
used this clean spectrum to interpolate linearly at a desired wave¬ 
length (the clean spectrum is highly sampled, so it is not necessary 
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Figure Al. In the left-hand panel, we show a CCD subsection with two night-sky emission lines. In the right-hand panel, we show two spectra for these lines: 
one from a single CCD row (red thick line), and other using all the rows of the CCD subsection (black thin line). Note that the spectrum from all the rows of 
the CCD subsection takes advantage of the C-distortion for getting a spectrum more sampled than the other from single CCD row. Bottom blue medium-width 
line shows the residuals in the CCD row of the red thick line after the subtraction of sky model constructed from all the rows of the CCD subsection. 



observed wavelength (A) 



observed wavelength (A) 


Figure A2. Top panel shows a CCD subsection; the two dashed lines indicate the region excluded for the sky modelling due to the presence of the target 
spectrum. Middle panel shows the spectrum build from the CCD subsection of the top panel taking advantage of the C-distortion. In bottom panel, we plot the 
sky model derived from the spectrum of the middle panel. 
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Figure A3. Top panel shows a raw CCD subsection. Middle panel shows the sky model for the CCD subsection. Bottom panel shows the CCD subsection 
after the sky subtraction. In each panel, dashed lines indicate the region near to the target excluded for building the sky model. Top and middle panels share 
the same colour bar (from 55 to 14409 counts), while another colour bar used for bottom panel (from 0.1 to 1000 counts). 


pyKelsame 



400 600 800 

spectral direction (rebinned pixel) 


Figure A4. Top/bottom panel shows a CCD subsection where the sky subtraction was performed before/after correcting the C-distortion and using the 
PYKelsame/SKYSUBM program. Both panels share the same colour bar (from 0.1 to 1000 counts). 
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a higher order interpolation). This sky model has been plotted in 
the bottom panel. 

Fig. I A3 1 illustrates process of the sky subtraction. The top 
panel is the raw CCD subsection used to build the sky model. The 
middle panel is the evaluation of the sky model in the pixels of 
the raw CCD. We can check how to the cosmic rays and the target 
spectrum have disappeared (as it was desired). The bottom panel 
shows the sky-subtracted CCD. Notice that the cosmic rays are 
still present in the CCD, but as the sky lines have been removed, 
it is easier to detect them (so afterwards it will be easier to remove 
them). 

Fig. I A4l compares the sky-subtraction method used in this pa¬ 
per with the traditional technique. Both of them have been per¬ 
formed in the environment of REDUCEME (using PYKELSAME and 
SKYSUBM programs, respectively). According to the above de¬ 
scriptions, the key difference between these methods is the order 
of subtracting sky and correcting C-distortion. For better compari¬ 
son, we have also corrected the C-distortion in the panel of Fig. lA4l 
which shows the sky subtraction with PYKELSAME: this allows 
that both panels display rebinned data. Checking the intensity of 
residuals after sky subtraction, Fig. I A4I confirms that the method 
explained in this appendix reports better results than the traditional 
technique. 


APPENDIX B: ROBUSTNESS OF THE VELOCITY 
DISPERSION MEASUREMENTS 

A number of tests indicate that the velocity dispersions and their 
errors are robust. 

We obtained velocity dispersions for different combinations 
of SSP templates. In particular, we repeated the PPXF fits: 

(i) Using the SSPs from the library with ages younger than age 
of the Universe at z = 1. 

(ii) Using the SSPs from the library with ages younger than age 
of the Universe at z = 1 and the option of regularization of PPXF. 

(iii) Using the whole SSP library. 

(iv) Using the whole SSP library, and masking H 7 and H<5 
lines. 

(v) Using the single SSP which alone best fits each spectrum. 

For all galaxies, the results from all the tests were in agreement 
with each other within their errors, being the mean standard devi¬ 
ation 9 km s” 1 . In the case (iii), we tested for the possibility that 
residual emission might be partially filling the H 7 and H <5 lines; 
because most of our spectra do not reach the [O II] 3727 A line, star 
formation could easily go unnoticed and affect the velocity disper¬ 
sion measurements. We therefore compared our results with those 
from fits in which these lines were masked out. 

We further tested out velocity dispersion results by carrying 
out the PPXF fits using a stellar library instead of SSPs. We used 
a library with 193 stars from the Indo-US Library of Coude Feed 
Stellar Spectra dValdes et al.li2004l) . The selection of stars was per¬ 
formed to cover a wide range of effective temperatures, metallici- 
ties and surface gravities. This library was chosen because its high 
spectral resolution (1.36 A FWHM) allows us to take advantage of 
the quality of our spectra (~2 A FWHM at rest frame). The stellar 
library has the disadvantage of being less physical motivated, but, 
on the other hand, it provides wider spectral coverage by reach¬ 
ing down to 3460 A, as compared to the PEGASE-HR SSP library, 
which cuts at 3900 A. These fits were carried out masking and with¬ 
out masking H 7 and H <5 lines. In all these cases, we obtain similar 


results: the mean of the differences is 0 and their standard deviation 
is 7 km s” 1 . 

The results from fitting stellar templates are very similar to 
those from SSPs for four of the galaxies. For the other two, the fits 
with stellar templates give lower, marginally inconsistent disper¬ 
sions. Specifically, for the galaxy ID 13018611, we obtained a = 
118±22kms” 1 using the stellar library and a = 151±15kms _1 
using SSPs (not masking the H 7 and H 5 lines in both cases). Sim¬ 
ilarly, for the galaxy ID 12028173, using the stellar library we ob¬ 
tained 180±8 and 214± 16 km s” 1 using SSPs (bootstrap errors). 
As discussed in Section [Til we believe the fits to stellar templates 
are more uncertain given the intrinsic width of the Ca II H and K 
lines for some stellar temperatures, and hence we choose the so¬ 
lutions given by SSP fits. Nevertheless, adopting the dispersions 
obtained with stellar templates would not change the conclusions 
from this work. 

Finally, we have compared our results with those from other 
authors. We have the galaxy ID 12024790 in common with MM11. 
That authors reports a e = 160± 10 km s” 1 , which is significantly 
lower than our value, a e = 261 ± 15 km s” 1 . Applying the mask¬ 
ing they use to exclude regions affected by telluric absorption, and 
cutting the redder part of our spectmm to their red cut-off we could 
not reproduce their result: we obtained a e = 258 ± 22 km s” 1 . 
In order to explore the origin of the difference, we repeated the 
measurement over the lMMl ll spectrum (albeit without relative-flux 
calibration) using PPXF and the above libraries (PEGASE-HR and 
Indo-US). Using the PEGASE-HR library, we performed the fit for 
a wide range (1-13) of degrees of the Legendre polynomial used to 
correct the continuum shape and found a strong dependence with 
this parameter obtaining values from 198 to 310 km s” 1 . Doing 
the same exercise with the Indo-US library, we did not find a clear 
dependence with the polynomial degree and the results cover a nar¬ 
rower range: 227-259 kms” 1 . We also performed these tests in 
our spectrum for this galaxy (using the same wavelength ranges 
and masking the same regions), and obtained the ranges 250-264 
and 239-277 km s” 1 for PEGASE-HR and Indo-US libraries, re¬ 
spectively. All these tests except one produce stable results similar 
to our measurement, while the unstable test corresponds to a sim¬ 
ilar situation to the published by iMMlH : their spectrum and SSP 
templates. Additional checks on our own spectra included a veri¬ 
fication of the instrumental dispersion and its variation from night 
to night: measurement of the width of the telluric lines in the re¬ 
gion of the kinematic signatures shows a very stable configuration 
(3.88 ± 0.06 .A FWHMs). We conclude that the origin of the dif¬ 
ference with ImmUI is unclear, and might be related to the use of 
different velocity dispersion algorithms; we have used our measure¬ 
ment in the paper. 

We share the galaxy ID 13018611 with the sample of 
iFernandez Lorenzo et al.l ( KOllh . For this galaxy, they report a e = 
196 ± 6 kms” 1 , while we find a e = 161 ± 16 kms” 1 . Their 
velocity dispersion measurement shows a dependence on the spec¬ 
tral range of the fit: 173 kms” 1 for the range 3829-4129 A, 
184 kms” 1 for 3775-4178 A and 195 kms” 1 for 3725^1227 
A (Fernandez Lorenzo, private communication). The first of these 
three ranges is the closest to the one used by us, and, for that fitting 
range, the two determinations agree within the errors. In any case, 
the differences found show that the true uncertainties in the veloc¬ 
ity dispersion determinations can be underestimated using the stan¬ 
dard techniques (Monte Carlo s imulations with white noise); this 
fact has been already proven by Ide Bruvne et al .i (2003 ). We have 
used our value in the paper for consistency, and emphasize that us¬ 
ing their value would not modify the conclusions of this work. 
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Figure Cl. pPXF fitting results for the targets in the slit 1 of the observations. On top panels black thin solid lines show the spectra, red thick solid lines 
represent PPXF fits, and blue dashed lines are lcr deviations from the PPXF fits. On bottom panels we have plotted the residuals with black dots. Units are 
arbitrary. Grey regions indicate those which were masked in the PPXF fits. We have marked some spectral features with green dashed vertical lines. 
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Figure C2. pPXF fitting results for the targets in the slit 2 of the observations. On top panels black thin solid lines show the spectra, red thick solid lines 
represent PPXF fits, and blue dashed lines are lcr deviations from the PPXF fits. On bottom panels we have plotted the residuals with black dots. Units are 
arbitrary. Grey regions indicate those which were masked in the PPXF fits. We have marked some spectral features with green dashed vertical lines. 
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Figure C3. PPXF fitting results for the targets in the slit 3 of the observations. On top panels black thin solid lines show the spectra, red thick solid lines 
represent PPXF fits, and blue dashed lines are la deviations from the PPXF fits. On bottom panels we have plotted the residuals with black dots. Units are 
arbitrary. Grey regions indicate those which were masked in the PPXF fits. We have marked some spectral features with green dashed vertical lines. 
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APPENDIX C: PPXF FITTING RESUFTS 

We show in Figs IC 1 ll('.3| the spectra obtained after the data reduc¬ 
tion explained in Section [U2l together with the spectral fit obtained 
from PPXF detailed in Section [33l 


APPENDIX D: STRUCTURAE PARAMETERS FOR THE 
Z r ^ J 1 SAMPLE 

Table iDll contains a compilation of structural parameters for the 
sample of massive ETGs at 2 ~ 1 analysed in this paper. 
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1 sample. (1) References: [0] this work; [11 [Martinez-Mansoetal. 


Table Dl. Compilation of structur al para met ers for the 2: __ _ 

1201 ll) ; [2][BeJli Ji hJewman&Ellis l2014|); [3llNewman et al. 1201 Oh (extracted from the compilation m ade by 
1201 


van de Sande et al. 

; [41 Ivan der Wei et alj 1200S ) and lBlakeslee et alj 120061) (extracted from the compilation made by Ivan de Sande et all2013i) . 


(2) Galaxy identifications as given in the references. (3) Spectroscopic redshifts. (4) Sersic indices. (5) Circular effective (half-light) 
radii. (6) Stellar masses (adapted to Salpeter IMF). (7) Velocity dispersions within the effective (half-light) radius. 


Reference 

ID 

z 

n 

r e 

M* 






(kpc) 

(10 11 M e ) 

(km s -1 ) 

a) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

0 

13004276 

0.97533 

5.1 

1.46 

2.0 

340 ± 16 

0 

12028173 

0.95764 

5.9 

1.42 

1.2 

228 ± 17 

0 

12024790 

0.96572 

4.4 

0.54 

0.9 

261 ± 15 

0 

12024422 

1.02874 

4.8 

1.57 

1.2 

239 ± 18 

0 

13017337 

0.97242 

4.3 

1.37 

1.8 

291 ±25 

0 

13018611 

1.07939 

4.7 

1.19 

3.0 

161 ± 16 

1 

12024321 

0.91592 

6.5 

2.46 

0.9 

162 ± 12 

1 

12019899 

0.93259 

6.3 

0.47 

1.1 

245 ± 18 

1 

12024453 

0.90564 

4.8 

1.60 

1.7 

184 ± 10 

2 

51106 

1.013 

5.3 

5.99 

3.3 

252 ± 37 

2 

28739 

1.029 

3.5 

1.98 

1.6 

238 ± 11 

2 

54891 

1.081 

3.2 

1.27 

1.0 

232 ± 37 

2 

31377 

1.085 

4.9 

4.88 

1.2 

133 ± 18 

2 

13393 

1.097 

3.5 

7.18 

2.6 

175 ± 21 

2 

16343 

1.098 

8.0 

1.95 

2.0 

290 ±8 

3 

El 

1.054 

4.0 

6.44 

2.1 

204 ± 22 

4 

761 

1.01 

4.0 

3.77 

6.2 

377 ± 40 

4 

1559 

0.94 

4.0 

1.66 

1.1 

179 ± 13 

4 

1706 

0.91 

4.0 

2.23 

2.7 

217 ± 13 

4 

1 

1.09 

4.0 

2.83 

2.7 

233 ± 17 

4 

2 

0.96 

4.0 

2.30 

3.5 

202 ± 10 

4 

3 

1.04 

4.0 

1.00 

1.2 

302 ± 34 

4 

4 

0.96 

4.0 

6.84 

7.8 

337 ± 19 

4 

13 

0.98 

4.0 

2.20 

3.2 

249 ± 11 

4 

14 

0.98 

4.0 

2.80 

1.4 

199 ± 24 

4 

18 

1.10 

4.0 

3.97 

5.7 

327 ± 37 

4 

20 

1.02 

4.0 

2.24 

2.6 

201 ± 17 
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